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Abstract 


Strong  Planning  and  Forecast  Horizons  for  a  Model  with 
Simultaneous  Price  and  Production  Decisions 

by 

G.L.  Thompson,  S.P.  Sethi,  and  J.T.  Teng 

In  this  paper,  we  have  solved  a  general  inventory  model  with 
simultaneous  price  and  production  decisions.  Both  linear  and 
non-linear  (strictly  convex)  production  cost  cases  are  treated. 
Upper  and  lower  bounds  are  imposed  on  state  as  well  as  control 
variables.  The  problem  is  solved  by  using  the  Lagrangian  form  of 
the  maximum  principle.  Strong  planning  and  strong  forecast 
horizons  are  obtained.  These  arise  when  the  state  variable 
reaches  its  upper  or  lower  bound.  The  existence  of  these  horizons 
permits  the  decomposition  of  the  whole  problem  into  a  set  of 
smaller  problems,  which  can  be  solved  separately,  and  their 
solutions  put  together  to  form  a  complete  solution  to  the  problem. 
Finally,  we  derive  a  forward  branch  and  bound  algorithm  to  solve 
the  problem,  the  algorithm  is  illustrated  with  a  simple  example. 
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Strong  planning  horizon 
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1.  Introduction 


Many  methods  have  been  developed  to  solve  the  problem  of 
production  scheduling  over  future  periods  of  time  in  order  to 
satisfy  the  given  amounts  required  of  a  certain  commodity  in  each 
of  these  periods.  Such  methods  are  based  on  a  wide  range  of 
production-inventory  models — from  elementary  deterministic  ones  to 
sophisticated  stochastic  models  with  time-varying  parameters. 
Historically,  applications  of  optimization  methods  to  production 
and  inventory  problems  date  back  at  least  to  the  classical  EOQ 

model  and  the  lot  size  formula.  The  EOQ  model  is  essentially  a 

static  model  in  the  sense  that  it  assumes  the  demand  to  be 

constant  and  only  a  stationary  solution  is  sought.  A  dynamic 
version  of  the  lot  size  model  in  which  the  assumption  of  a 
steady-state  demand  rate  was  dropped  was  analyzed  by  Wagner  and 
Whitin  [13].  The  solution  methodology  used  there  was  dynamic 
programming . 

In  some  dynamic  problems  it  is  possible  to  show  that  the 

optimal  decisions  made  during  an  initial  positive  time  interval 
are  either  partially  or  wholly  independent  of  the  data  from  some 
future  time  onwards.  In  such  cases,  a  forecast  of  the  future  data 
needs  to  be  made  only  as  far  as  that  time  to  make  optimal 
decisions  in  the  initial  time  interval.  The  intial  time  interval 
is  called  the  planning  horizon  and  the  time  up  to  which  data  is 
required  to  make  the  optimal  decisions  during  the  planning  horizon 
is  called  the  forecast  horizon.  Whenever  they  exist,  these 
horizons  naturally  decompose  the  problem  into  a  series  of  smaller 
problems . 
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If  the  optimal  decisions  during  the  planning  hcrison  are 
completely  independent  of  the  data  beyond  the  forecast  fcrizcn, 
then  the  former  is  called  a  strong  planning  horizon  and  the  latter 
is  called  a  strong  forecast  horizon.  If,  on  the  other  hand,  some 
mild  restrictions  on  the  data  after  the  forecast  horizon  are 
required  in  order  that  to  keep  the  optimal  decisions  during  the 
planning  horizon  unaffected,  then  it  is  called  a  weak  forecast 
horizon ,  and  the  corresponding  planning  horizon  is  called  week 
planning  horizon.  For  example,  Modigliani  and  Hchn  [5]  shew  that 
a  planning  horizon  [0,t  ]  in  their  production  planning  model,  the 
corresponding  strong  forecast  horizon  is  always  the  overall 

4 

horizon  [0,T],  but  the  information  required  on  [t  ,T]  is  just  the 
accumulated  demand  rather  than  the  complete  demand  schedule  for 
each  instant  of  time  in  that  interval.  In  the  case,  the  weak 
planning  horizon  and  the  weak  forecast  horizon  are  identical.  Of 
course,  a  weak  planning  horizon  does  not  reduce  the  interval  over 
which  forecasts  are  needed,  but  reduces  the  informational 
requirement  within  this  interval. 

Earlier  planning  horizon  results  for  production  planning, 
obtained  by  using  dynamic  programming  and  variational  arguments, 
have  been  given  by  Kunreuther  and  Morton  C 2 , 3 D *  Modigliani  and 
Hchn  [5],  and  Wagner  and  Whitin  [13].  Other  horizon  results, 
derived  by  using  optimal  control  theory,  were  given  by  Kleindcrfer 
and  Lieber  [1],  Lieber  [4],  Morten  [6],  Pekelman  [7,8],  Teng, 
Thompson,  and  Sethi  [10],  and  Vanthienen  [12]. 

In  this  paper,  the  problem  of  determining  simultaneously  the 
price  and  production  schedule  of  a  firm  over  a  finite  horizon  T  is 


solved  for  a  given  time  dependent  demand  function.  Beth  linear 
and  non-linear  (strictly  convex)  production  cost  cases  are 
treated.  Upper  and  lower  bounds  are  imposed  on  state  as  well  a3 
control  variables.  The  problem  is  solved  by  using  the  Lagrangian 
form  of  the  maximum  principle.  Strong  plai.ning  and  strong 
forecast  horizons  are  obtained.  These  arise  when  the  state 
variable  reaches  its  upper  or  lower  bound.  The  existence  of  these 
horizons  permits  the  decomposition  of  the  whole  problem  into  a  set 
of  smaller  problems,  which  can  be  solved  separately,  and  their 
solutions  put  together  to  form  a  complete  solution  to  the  problem. 
Finally,  we  derive  a  forward  branch  and  bound  algorithm  to  solve 
the  linear  case.  The  forward  algorithm  for  the  convex  case  is 
similar. 

2.  The  Linear  Production  Cost  Case 

» 

In  this  section,  we  will  find  the  optimal  price  p  (t)  and 

ft 

production  rate  u  (t)  of  a  monopolist  who  faces  a  differentiable 
demand  curve  on  the  interval  [O.T].  Suppose  the  demand  curve  is 
D(t)  r  a(t)-bp(t)  and  the  linear  production  cost  is  cu(t),  the 
linear  inventory  cost  is  hl(t),  the  maximum  production  rate  is  u, 
and  the  warehouse  constraint  is  I(t)  <  W.  The  resulting  problem 
is  that  of  maximizing  the  following  expression: 

Jc  =  /q(p( t) C  a ( t)-bp( t) ]-hI ( t)-cu( t) }dt  (1) 


subject  to 
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l\l  i(t)  =  u(t)-[a(t)-bp(t) ],  1(0)  =  IQ  <  W  (2) 

[p,]  I(t)  >  I(t)  >  0,  [p2]I(t)  <  W,  [ u 1 ] u( t)  >  0, 

(3) 

Cu2]  u(t)  <  u,  tn^  p(t)  >  0,  t n2 ] pC t >  <  a(t)/b, 

where  h,c,b,u  and  W  are  positive  constants  and  c<a(t)/b  for  all  t. 
Note  that  a(t)/b  represents  the  maximum  price,  which  is  required 
to  exceed  the  production  cost  (otherwise,  no  production  will  take 
place).  The  function  x(t)  is  the  adjoint  variable  of  (2),  and 
Pl(t),  p2(t),  Ml(t),  U2(t),  n.|(t)  and  are  the  Lagrange 

variables  of  the  corresponding  constraints.  A  dot  above  a 
variable  denotes  the  first  derivative  with  respect  to  time. 

This  model  is  similar  to  the  model  due  to  Pekelman  [7]  except 
that  the  production  cost  here  is  linear,  and  production  and 
warehouse  upper  bound  constraints  have  been  added.  In  Section  3 
we  consider  the  same  model  with  convex  production  cost,  which  is  a 
generalization  of  Pekelman's. 

2.1  The  Necessary  Conditions  for  an  Optimal  Solution 
The  Hamiltonian  function  of  Problem  (1)  is 

H=pa-p^b-hl-cu-* (u-a+bp)  (4) 

which  is  linear  in  u  and  quadratic  in  p  sc  that  the  optimal 
production  rate  u*(t)  is  bang-bang  and  the  optimal  price  p#(t)  is 
a  sat  function.  For  this  problem  the  Lagrangian  is 


t.»Ht»1(u.a*bp)+,2(a.bp.a)tll)11„2(-_u)+11^ii^(>/b_ii)  (5) 

Th«  following  necessary  conditions  hold,  see  »],  for  an  optimal 
solution 


AL 


u  ^-c-^+vvw 


or 


* 

u  .< 


undefined 
l  u 


if  X+Pi_p2<c 

if  x+p^p^c 

if  x+Pi_p2>c; 


(6) 


AL 


ap  =0sa-2bp+(i1+P1-p2)b+n1-n  cr 


r  0 


* 

p  M 


if  x+Pi_p2<_a/b 


(a/b+x+P 1-P2)/2  if  -a/b<^+P 1-02<a/b 

(.a/b 


if  a/b< i+P i -P2 ; 


(7) 


the  adjoint  equation  satisfies 


x=-3L/3I=h 


(8) 


and  the  transver'sality  conditions  are 


x (T) >0  and  x(T)I(T)=0; 


(9) 


i 


the  complementarity  and  nonnegativity  conditions  are 


p-|I»p1I,  p2(W-I  ) ,  Pgl  »uiu,m2(u-u)  ,  n2(a/b-p)  sO 


(10) 


L  and  p1.p2*ui.W2*ni  .n2>°; 

■the  variable  I(t)  is  continuous,  x  is  continuous  except 
possibly  at  an  entry  time  or  an  exit  time  to  the  boundary 

conditions  I(t)s0  or  I(t)=W,  x+p^-p2  xs  continuous  everywhere, 

•  • 

land  moreover  P1^0  and  p2<0.  (11) 

Note  that  t.  is  an  initial  or  entry  time  to  I(t)  =  W  if  I(t.  )  < 

1W  ■—  ■  — 1  •  - - *-  1W 

W,  and  I  ( tiw)  =1  ( tX)  =W;  tfw  is  a  final  or  exit  time  to  I(t)=W  if 
I  (tf~)=I  (tfw)=W  and  I(tf^)<W.  The  definition  of  tiQ  and  tf0  are 
similar  and  are  obtained  by  changing  W  to  0. 

2.2  Optimal  Policies  for  Three  Possible  Cases 

There  are  only  three  different  possible  cases  for  the  values 
of  I(t).  We  shall  discuss  the  optimal  policies  in  cases  1,2, and  3 
below . 

Case  1  0<I(t)<W 

This  implies  p^sP2=0.  We  have  two  different  regimes  1A  and 
IB  in  this  case. 

Regime  1A  x(t)<c.  From  (6)  and  (7)  we  have  the  following 
optimal  policies; 

u  =0,  p  =max l0 ; (a/b+x )/2 } ,  I<0 


(12) 
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Regime  IB  x(t)>c. 

u*=u  and  p*=min  (a/b;  (a/b  +  x)/2*  (13) 

Case  2  l(t)=0  and  l(t)=0. 

This  implies  p2=0,  u  =a-bp  ,  and  ( a/b+c)  /2<p* <a/b.  We  also 
have  regimes  2A  and  2B  in  this  case. 


Regime  2A  u>(a-bc)/2.  An  application  cf  the  maximum 
principle  gives  the  following  optimal  policies; 

u*= (a-bc) /2 ,  p*= (a/b+c) /2 ,  x+p1=c  and 

(14) 

",  *  "2  =  "l  ’  "2  =  ° 

Regime  2B  u<(a-bc)/2  and  a<bh.  The  maximum  principle  give 
the  following  optimal  policies; 


u*=u,  p  = ( a-u ) /b ,  x+p ^ = ( a-2u ) /b , 
y 1 =n 1 =^2=0,  and  P2= ( a-2u-bc )/b>0 


(15) 


•  • 

Note  that  if  both  u<(a-bc)/2  and  a>bh  hold  then  p1=a/b-h>0,  which 
leads  to  a  contradiction  to  (11). 


Case  3 


I(t)sW  and  l(t)=0. 
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These  conditions  imply  ^=0,  u  =  a-bp  , 
We  have  only  one  regime  3A  in  this  case. 


# 

and  ( a/b+c) /2<p  <a/b. 


Regime  3A  u<(a-bc)/2  and  a>bh.  The  maximum  principle  yields; 


u*=u ,  p*=  (  a— u)/b  ,  X-P2=  ( a-2u)/b  ,  ^  ^  =  1  =  n2  =  0, 


and  »25 ( a-2u-bc )/ b>0 


(15) 


*  * 

Note  that  if  u>(a~bc)/2  then  u  =(a-bc)/2,  p  =(a/b+c)/2,  and 

•  •  •  • 
i-P2=c.  This  implies  X-P2=0  and  contradicts  the  fact  that  A-P?>h. 

•  „  •  • 

Similarly,  if  u<(a-bc)/2  and  a<bh  then  x-P2=(a-2u)/b,  X-P2=a/bsh 
and  this  again  leads  to  a  contradiction.  Therefore  Regime  3A  is 
the  only  case  to  be  considered  here. 


A  simple  economic  interpretation  of  Cases  2  and  3  is  as 
follows.  It  is  easy  to  prove  (proof  omitted)  that  1=0  implies 
that  (a-bc)/2  is  the  best  production  rate  without  considering  the 
production  capacity.  Therefore,  we  can  call  (a-bc)/2  an  ideal 
perfect  production  rate.  In  Regime  2A,  u>(a-bc)/2  means  that  we 
can  3et  the  production  rate  to  be  the  ideal  perfect  production 
rate  without  violating  the  production  capacity.  Thus,  we  don't 
need  to  build  up  inventory  for  the  future,  i.e.,  Keeping  the 
warehouse  empty  is  cur  optimal  strategy.  In  Regime  2B,  u<(a-bc)/2 
implies  that  the  production  rate  couldn't  reach  to  the  ideal 
perfect  production  rate  because  of  a  lower  production  capicity. 
However,  a<bh  hints  that  the  rate  of  price  increase,  a/b,  is  net 
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greater  than  the  marginal  inventory  cost,  h.  Therefore,  it  does 
net  pay  to  build  up  inventory.  Similarly,  in  Regime  3A  the 
production  capacity  is  less  than  the  ideal  perfect  production  rate 
and  the  rate  of  price  increase  is  greater  than  or  equal  tc  the 
marginal  inventory  cost.  Thus,  it  pays  to  build  up  inventory. 

2.3  Theoretical  Results 

Let  us  investigate  the  behavior  of  I(t).  First,  we  shall 

establish  a  sufficiently  large  initial  inventory  IQ  such  that 

_  * 

Iq>Iq  implies  u  =0  for  all  t.  From  Regime  1A,  and  Equation  (9)  we 
knew 

I0  =  'J[a-bmax{0;  [a/b+h( t-T ) ] /2 >  ]dt  (17) 

The  inventory  IQ  is  just  sufficient  tc  meet  demand  when  price  is 
determined  by  either  of  the  first  two  rules  in  (7). 

Theorem  1  Iq>T0<=>I (T) =Iq-T0>0,  and  *(T)=0 

<=>u  =0  for  all  t,  and  p„ (T ) =a/ (2b) . 

#  „ 

Proof.  We  prove  the  second  implication  only.  Assume  u  =0 
and  p  (T)  =  a/(2b).  This  can  occur  only  in  Regime  1A.  Thus, 
0<I(t)<W,  A(t)=h(t-T)  for  all  t,  and  I  (T )  si  -I  >0.  The 

verification  of  the  converse  is  trivial. 


Corollary  1  Iq<  Iq<=  >1  (T )  =0  and  X(T)>0 
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# 

Theorem  2  We  can  conclude  that  I<0  under  each  of  the  following 
cases : 

ro-*o 

IQ<I0  and  u>(a-bc)/2  for  all  t 
Iq<Iq;  whenever  u<(a-bc)/2  we  have  a<bh. 

Proof .  Case  (i)  follows  from  Theorem  1.  We  now  prove  cases 
(ii)  and  (iii).  Suppose  not,  i.e.,  I(t)>0  in  (t,t+c);  see  Figure 

1.  This  can  happen  only  in  Regime  IB.  Thus,  for  all  te(t+e/2, 

*  * 

r  +  e)  we  have  x>c,  u  =u  and  p  =min[a/b;  (a/b+x)/2].  Hence 

.  _  »  fix  if  X>a/b  (p*=a/b) 

I  =  u-( a-bp  )  = j  * 

lu-(a-b*)/2  if  Ua/b  (p  =(a/b+X)/2) 

In  Case  ( ii)  ,  the  conditions,  u>(a-bc)/2  with  x(T  +  e)>c,  imply 
u-(a-bi )/2>u-(a-bc)/2>0  on  (t+e/2,  t+e).  Thus  I(t+e)>0.  By 

repeating  this  argument  for  the  intervale  (t+e,  t+2e),  (t+2e, 

T+3e),  .....  we  see  that  I(T)>0.  This  ccntraditics  Corollary  1 
which  says  l(t)=0  in  Case  (ii). 

In  Case  ( iii ) ,  1=0  or  (bh-a)/2>0.  Thus,  i(T+e)>o,  which,  as 

above,  contradict  I(T)=0  for  Case  (iii). 


(i) 

(ii) 

(iii) 


T  X+t/j.  X +£ 


Figure  1 

Corollary  2  Under  the  cases  specified  in  Theorem  2,  the  optimal 
path  for  I(t)  is  in  one  of  the  following  three  cases: 

(a)  I ( t) >0  for  all  t  (This  case  can  happen  if  and  only  if  I0>TQ) 

(b)  l(t)=0  for  all  t  (This  case  can  happen  if  and  only  if  IQ=0) 

(c)  l(t)>0  on  C 0 , t ^ 3  and  l(t)=0  on  [t^tT]  (This  case  can  occur  if 

and  only  if  0<In<7„) 
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We  new  discuss  hew  in  Corollary  2(c)  is  determined.  For 
convenience,  let  us  define  the  function  <|»  ( t)  =  (  x+p  ^-Pg)  ( t)  which  is 
similar  to  the  function  in  Pekelman  [7].  In  Regimes  2A  and  2B, 
we  have 


{c  if  u>(a-bc)/2 

(18) 

(a-2u)/b  if  u<(a-be)/2  and  a<bh 

or  <l>Q=max(c;  (a-2u)/b).  Intervals  such  that  u>(a-bc)/2  or 
tu<(a-bc)/2  and  a<bh]  are  defined  to  be  empty  warehouse  transitive 
intervals  on  which  the  boundary  condition  l(t)=0  holds. 
Similarly,  by  using  the  results  in  (16)  for  Regime  3A  we  have 

( a-2u)/b  (19) 

w 

Intervals  such  that  u<(a-bc)/2  and  a>bh  are  defined  to  be  full 
warehouse  trasitive  intervals  on  which  the  bcudnary  condition 
I(t)=W  holds. 

After  defining  the  function  1<q,  we  can  use  a  binary  search 
program  to  find  the  values  of  x(0)  and  t1  in  Corollary  2(c).  By 
Corollary  1  we  have  iQsX (0 ) >-hTsXQ .  If  condition  (ii)  of  Theorem 
2  occurs  then  *0(t1)=c  sc  that  xg<°.  If  condition  (iii)  of 
Theorem  2  happens  the  xQ<a(0)/b;  otherwise  X(t)>a(t)/b  (because 
x=h>a/b),  p(t)  =  a(t)/b  for  all  t  so  that  I(T)>IQ  and  leads  to  a 
contradiction.  Thus,  Xg<a(0)/b=XQ  holds  in  any  of  the  three 
conditions  of  Theorem  2.  To  find  the  value  of  t^,  we  can  first 
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guess  Xq^  =  (Xq+Xq)/2.  Check  Whether  xg^  satisfies  the  following 
3  equations  or  not. 


X(t1)SX0+ht1S*g(t1) 

(20) 

X(tc)  =  Vhtc=c 

(21) 

v  -i 

Ig  s-^g  ^a-bmaxCO;  (x0+ht+a/b)/2  ]Jdt-^t  (u- [a-b(  xQ+ht)  ]/2 >dt  (22) 

c 

If  yes,  then  we  are  dene.  If  the  right-hand  side  of  (22)  is 
large.r  than  I  ^  then  let  the  new  guessing  value  of  x  Q  is 
Xg^ ( Xg ^ +Xq) /2  because  I  is  an  increasing  function  with  respect 
to  xQ.  Otherwise,  let  xq2^  =  (Xq1 ^+Xq)  •  Repeating  this  binary 
search  procedure  we  can  get  the  values  of  t^  and  x^.  Next,  we 
shall  explore  the  case  in  which  a>bh,  i.e.,  the  case  in  which  it 
pays  to  store  inventory. 

Lemma  1 .  If  a  is  continuous  and  non-decreasing,  then  the  optimal 

value  of  +  is  also  continuous  and  non-decreasing. 

Proof.  In  Regimes  1A  and  IB,  * =x  is  an  increasing  continuous 

function.  In  Regimes  2A  and  2B,  =*0=max  (c ;  ( a-2u)/b)  is  a 

non-dec  •'easing  continuous  function.  Similarly,  in  Regime  3A 

*  =  *  =(a-2u)/b  has  the  same  properties, 
w 

Corollary  3.  If  a  is  continuous  and  non-decreasing  then  the 

ft  ft 

optimal  controls  u  and  p  are  both  non-decreasing. 
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Prccf .  It  immediately  follows  from  Lemma  1  and  (6),  and  (7). 

*  o  %  — 

Corollary  4.  If  a  is  continuous  and  non-decreasing  then  u  (t  )=u 
implies  u*(t)=u  for  all  t>t°. 

Proof .  Follows  immediately  from  Corollary  3  and  the 
constraint  u<u. 

Lemma  2 .  If  a>bh  then  l(t)=0  on  any  interval  (t^t1)  will  imply 
u*s(a-bc)/2<u  on  (t^t1). 

Proof.  I(t)=0  on  (t^t1)  implies  this  case  can  happen  only 
in  Regime  2A.  Thus,  u>(a-bc)/2  and  u*=  ( a-be)  /2<u  on  (t^t1). 

From  Corollary  4  and  Lemma  2,  we  can  easily  get  the  following 
result . 

Corollary  5.  If  a>bh  and  u  (t  )  =  u  then  l(t)>0  almost  everywhere 
on  [t°,T]. 

Theorem  3.  If  a>bh  then  there  exists  an  optimal  solution  such 
that  l(t)>0  almost  everywhere  on  [0,T]. 

ft  _ 

Proof.  To  prove  this  theorem,  we  change  (6)  sc  that  u  =u 

when  i+P^-p2^c>  which  comes  by  putting  together  the  last  two  parts 

* 

of  (6).  If  u  (t)=0  for  all  t  then  we  are  in  Regime  1A  so  that 

I ( t) >0  almost  everywhere  on  [0,T].  Otherwise,  there  exists 

*  0  -  0 

u  (t  )=u  for  some  t  *[0,T].  Let  us  define 
tQsinf {tiu*(t)=u> . 
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If  tQ=0  then  the  theorem  is  true  by  Corollary  5.  If  fco*°  and 
there  exists  an  open  interval,  say  (tg.t^)  with  t^tg*  such  that 

l(t)=0  on  C  t2 1 1 3 ) »  then  *g(t)>c  and  cn  (t2»fc3^’  This 

contradicts  the  fact  that  u  (t)=0  on  (t ^ ,  t^).  Thus,  l(t)>0 

almost  everywhere  cn  [0,tQ].  By  Corollary  5,  we  also  have  l(t)>0 
almost  everywhere  on  [tQ,T].  This  completes  the  proof. 

Theorem  4.  If  l(t)>0  almost  everywhere  cn  [0,T],  then  the  optimal 
control  path  for  u*(t)  satisfies  one  of  the  following  three  cases: 

(I)  u*=0  for  all  t 

(II)  u*=u  for  all  t 

«  #  _ 

(III)  u  =0  for  t<tQ  and  u  =u  for  t>tQ 

Proof.  If  I( t) >0  almost  everywhere  on  [0,T]  then  case  2 
cannot  occur  so  that  there  are  no  jumps  in  *(t),  except  in  Regime 
3A.  Regime  3A  can  occur  only  following  Regime  IB;  also  1A  cannot 
follow  Regime  IB.  Therefore,  we  have  the  following  three  possible 
optimal  policies: 

(1)  Regime  1A,  fc.  all  t.  (This  implies  (I)); 

(2)  Regime  IB  at  first,  then  Regime  3A  and  IB  maybe 

alternating  several  times,  at  time  T  the  case  should  be  Regime  IB. 

(This  implies  (II)); 

(3)  Regime  1A  at  first,  then  Regime  IB,  after  this  Regime  3A 
and  IB  maybe  alternating  several  times,  at  time  T  the  case  should 
be  Regime  IB.  (This  implies  (III)). 

This  completes  the  proof  of  Theorem  4. 

Here,  we  shall  also  discuss  the  situation  in  which  the 
exogenous  variable  a(t)  can  first  go  up  and  then  down  or  first 
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down  then  up.  In  fact,  using  the  above  results  we  can  easily 
obtain  the  following  theorems. 

Theorem  5 .  Assume  t  is  chosen  so  that  a>bh  for  t<t  and  a<bh  for 
t>t;  i.e.,  it  pays  to  store  production  before  but  not  after  t. 
Also  assume  neither  condition  (i)  nor  condition  (ii)  of  Theorem  2 
holds.  Then  there  exists  an  optimal  solution  such  that  I(t) 
satisfies  one  of  the  following  two  cases: 

(i)  I(t)>0  for  all  t«[OfT) 

(ii)  there  exists  t2>t  such  that  l(t)>0  on  (0,t2)  and  l(t)=0 
cn  C 1 2 » T ] . 

Theorem  6.  Assume  t  is  chosen  so  that  a<bh  for  t<t  and  a>bh  for 
t>t.  Also  assume  neither  condition  (i)  nor  condition  (ii)  of 
Theorem  2  holds.  Then  there  exists  an  optimal  solution  such  that 
I ( t )  satisfies  one  of  the  following  three  cases: 

(1)  l(t)>0  for  all  t e [ 0 , T) 

(2)  there  exists  t3<t  such  taht  I(t)=0  [0 , t 3 ]  and  l(t)>0  cn 
( 1 3 » T )  • 

(3)  there  exists  t1+<t^<t  such  that  l(t)>0  on  [0,tjj),  l(t)=0 
on  Ct^,tg],  and  l(t)>0  cn  (t^,T). 

In  general,  the  interval  [0,T]  will  contain  many  subintervals 
which  u>(a-bc)/2  or  a<bh  or  a>bh.  By  repeatedly  applying  the 
results  of  Corollary  2  and  Theorems  3>  5,  and  6  we  can  construct 
the  solution  by  piecing  together  different  solutions  obtained  from 
application  of  these  theorems. 
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2.4  Planning  Horizon  Theorem 

The  optimal  decisions  (price  and  production  rate)  during  the 
planning  period  [0,t  ]  can  be  completely  independent  of  the  data 

ft  ft  ft  ft  ft 

beyond  the  forecast  horizon  [0,t  ]  with  t  <t  .  That  is,  no 

ft  ft 

information  after  t  is  required  for  making  optimal  decisions  on 

ft  ft 

[0,t  ].  We  therefore  call  t  to  be  a  strong  planning  horizon  and 

*ft 

t  to  be  a  strong  forecast  horizon. 

In  this  section,  we  will  explore  the  strong  planning  and 
forecast  horizons  for  the  problem. 

ft 

Lemma  3  Assume  we  have  the  optimal  trajectory  such  that  t  and 

ft  ft 

t  are  the  two  consecutive  entry  times  at  which  inventory  hits 

ft  ft  ft 

the  boundary  ccntraint.  Furthermore,  let  I ( t  )=0  and  I(t  )=W,  or 

I(t*)=W  and  l(t**)  =  0.  Suppose  that  t*<t**.  Then  t*  is  a  strong 

«« 

horizon,  and  t  is  a  strong  forecast  horizon. 

»  . 

Proof .  We  may  assume,  without  less  of  generality,  that  t  is 
the  first  time  to  enter  the  boundary  condition  l(t)=0  or  W.  Let 

the  optimal  value  of  1(0)  be  1Q  to  the  original  problem.  Then  the 

*  ft 

values  of  t  and  XQ  are  decided  by  the  demand  in  [0,t  ],  the 

production  capacity,  and  the  initial  inventory  as  follows. 

r'“T . “ 

*•  0(or  W)-IQx/£  Idt . 
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* 

If  we  can  prove  that  i(t)=iQ+ht  on  (0,t  )  is  still  our  optimal 

solution  to  any  problem  having  the  same  information  as  the 

»*  »• 

original  problem  on  [0,t  ]  and  regardless  of  its  value  after  t  , 

then  we  are  done.  Here,  we  will  prove  the  case  in  which  I(t*)=W 
and  I(t  )  =  0,  see  Figu.  e  2.  In  fact,  using  an  analogous  argument 
we  can  prove  the  other  case  in  which  l(t*)=0  and  I(t**)=W. 


If  the  new  optimal  value  of  1(0)  were  iQ  with  x0>xQ,  then  the  new 

* 

inventory  I  should  be  greater  than  the  original  I  on  [0,t  J 
because  the  higher  value  of  *  implies  the  higher  value  of  I  by 

ft 

using  (2),  (6),  and  (7).  Thus,  there  exists  a  te(0,t  )  such  that 
I(t)sW.  Let  t°  be  the  earliest  time  such  that  I(t)=W.  Then  we 
have  l(t)>0  for  all  te(t°,t*).  Otherwise,  let  t1  be  the  first 
time  such  that  Kt^sO  and  t^eCt^t*),  see  Figure  3  for 
illustration.  This  implies  0<I<W  and  ^hs*  because  Regimes  2A  and 

A  i 

2B  cannot  occur  on  (t  ,t  )  for  I. 
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Figure  3 


Since  T(  t°)  ='xo+ht0>xo+ht°  =  4'  ( t 0 )  ,  H>h  =  i  cn  [t0,^]  we  have  I>I  cn 
C  t 0 , t 1 1  and 

I(t1)-I(t°)=^tQidt<^toIdt=I(t1)-I(t0)=-W, 
t  t 

cr  Kt1  )=I  (t°)-W<0 ,  which  leads  tc  a  contradiction.  Thus,  l(t)>0 

for  all  te(t°.t*).  Again,  0<I<W  on  (t°,t*)  implies  T<>h=i  and  ¥>i 
,  0  * 

cn  (t  ,t  ).  Then  we  have 

nt  )-W=/t  Xdt>/tnidt=W-I(t°)>0, 
tu  "  tu 

which  again  leads  tc  a  contradiction  tc  the  fact  that  I(t)<W  for 

/V 

all  t.  Therefore,  we  have  shown  that  xg>X0  cannot  be  an  optimal 

solution  tc  the  new  problem.  Next,  if  the  new  optimal  value  of 

# 

x(0)  were  with  ^g<x0.  then  the  new  inventory  ^<1  cn  [0,t  ]  sc 
that  £(t  )<W,  see  Figure  4.  Since  ♦  sfw  on  (t  ,t^)  where  t^  is  the 
final  cr  exit  time  of  I(t)=W  we  know  a>bh,  ,j<h<4',  .£<♦,  £<i,  and 
£<1  cn  (t*,tf).  Similarly,  we  can  get  jK*,  <i,<X ,  and  £<W  on 
[tf,t**].  This  implies 
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or  i 


I(t**)-I(tf)<I(t**)-I(tf)=-W, 


(t 


i* 


)<0  and  leads  to  a  contradiction  again 


I 
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I 


Theorem  7.  (Strong  Planning  Horizon  Theorem) 

«  f 

Let  t  and  t  be  the  entry  and  exit  times  to  I(t)=W 

respectively,  i.e.,  I(t)=W  for  te[t  ,tf],  and  I(t  “)  and  I(tf+)<W. 

»« 

Suppose  t  is  the  next  entry  time  to  I(t)sO.  Then  all  t  before 

f  »* 

t  are  strong  planning  horizons  and  all  t  after  t  are  strong 

f 

forecast  horizons,  i.e.  t  is  a  maximal  strong  planning  horizon 

ft  * 

and  t  is  a  mimimal  strong  forecast  horizon.  The  theorem  also 
holds  when  0  and  W  are  interchanged. 

f 

Proof.  If  t  is  not  a  strong  planning  horizon  then  there 

0  *  f 

exists  a  new  optimal  solution  such  that  t  e(t  ,t  )  is  the  new  exit 

from  I ( t) =W  boundary.  The  reasons  are  as  follows.  If  t°<t  then 

« 

it  contradicts  to  Lemma  3  that  t  is  a  strong  planning  horizon. 
On  the  other  hand,  if  t°>tf  then  tf  is  a  strong  planning  horizon. 
Using  the  same  arguments  as  in  Lemma  3  we  can  prove  that  the  new 
exit  time  t°  will  lead  to  a  contradiction. 

2.5  A  Forward  Branch  and  Bound  Algorithm 

The  solution  to  (1)  is  obvious  if  I0>IQ.  Therefore,  we  may 
assume,  without  loss  of  generality,  that  I0<IQ.  In  order  to  find 
an  upper  bound  to  the  optimal  objective  value  of  (1)  we  solve  the 
following  problem: 

max  Jt=/^(p(t) [a(t)-bp(t) ]-cu(t) >dt,  (23) 


subject  to 


i(t)=u(t)-[a(t)-bp( t) ],  I(t)=It 


(24) 


I 
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It  is  easy  to  shew  that  the  optimal  solutions  are  given  by: 

P* ( t) = ( a+bc) / (2b) ,  D* ( t ) = ( a-bc) /2 

and 

Jt=/^(a2-b2c2)/(4b)dt-c[/J(a-bc)/2dt-It]  (25) 

For  convenience,  we  let 

Q 

J[a,8]smax/  [p( a-bp) -hl-cu]dt ,  subject  to  constraints  (2)  and  (3)» 

It  is  clear  that  Jt  in  (25)  is  an  upper  bound  of  J[t,T]  because  we 
do  net  have  inventory  cost  and  production  capacity  constraint  in 
problem  (23).  Next,  we  shall  say  that  the  vertex  v.  is  f athemed  if 

« J 

and  only  if  no  further  exploration  from  this  vertex  can  be 
profitable.  Otherwise,  we  shall  say  that  v^  is  unf athemed  or 
alive .  We  are  new  in  a  position  to  present  the  algorithm. 

Forward  Branch  and  Bound  Algorithm 
Step  0  (Initalizaticn)  Begin  at  the  live  vertex  v^,  where  J=0.  Go 
to  Step  1 . 

Step  1  (Branching)  Assume  that  t^  is  the  first  entry  time  to  the 
constraint  I(t)sW  or  l(t)=0.  Solve  the  following  two  cases. 

Case  1 . 1  Suppose  that  tj  is  the  entry  time  of  I=W.  Then 
should  satisfy  the  following  constraints: 

x(t^)=[a(t^) )-2u]/b  if  u<(a-bc)/2  and  a>bh  at  t=tj  (26) 
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Mt)  =  Mt“)+h(t-t.) 

for  all  t<t^ 

(27) 

tf-I0>/J[u-(a-bp)]dt>-I0 

for  all  t<t^ 

(28) 

W-Iqs/q  [u-(a-bp)]dt 

(29) 

If  there  exists  one  or  more 

solutions , 

then 

we  keep 

each 

such 

solution  as  successor  vertices  of  Vq 

and 

go  to 

Case 

1 .2. 

Otherwise,  there  are  no  solutions  for  this  branch,  terminate  this 
branch  and  go  to  Case  1.2. 

Case  1  .2  Suppose  that  ti  is  the  entry  of  1=0.  The  ti 
satisifies  the  same  ccntraints  as  in  Case  1.1,  except  that  (26) 
and  (29)  must  be  replaced  by  (30)  and  (31),  respectively. 


c 


if  u>( a-bc) /2  at  tstt 


Ut') 


l[ a( t^ )-2u]/b  if  u<(a-bc)/2  and  a<bh  at  t=t^ 


(30) 


-Ig^Q^Cu-ta-bp)  ]dt  (31) 
Again,  we  keep  all  solutions  as  successor  vertices  of  vQ,  if  any, 
and  go  to  Step  2. 

Step  2  (Update  Bound)  Check  each  new  live  vertex  Vj.  If  t^^=T 
and  J^<J,  where  t^  is  the  last  time  of  Vj  at  which  I  reaches 
to  W  or  0,  and  J  ^  is  the  objective  value  of  v^,  then  the  vertex 
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Gc  tc 


v.  is  fathomed.  If  t(:i)=T  and  J(^>J,  then  let 
Step  3- 

Step  3  (Fathoming  by  Bound)  Check  each  new  live  vertex  v^.  If 
J[0, t^ j ^ ]+Jt^ j,  then  the  vertex  v^  is  fathomed.  Go  to  Step  4. 

Step  4  (Branching)  If  no  live  vertices  exist,  go  tc  Step  7; 
otherwise,  select  a  live  vertex  v..  If  I(t^^)=W,  then  go  tc  Step 

J 

5;  otherwise,  go  to  Step  6. 

Step  5  Assume  that  t^w  is  the  exit  time  from  I(t)=W  and  t^  is  the 
next  entry  time  of  1=0  or  I=W.  As  in  Step  1,  we  have  two  cases. 

Case  5.)  Suppose  t.  is  the  entry  time  tc  1=0  and  t-  <t.. 
——————  i  i  w  1 

Solving  the  following  ccntraints,  we  can  find  the  values  cf  fcfw 
and  t^. 


X(tfw*>  = 

=[a(tfw)-2u]/b 

if  u<(a-bc)/2  and  a>bh  at  t  =  t^w 

(32) 

|c 

if  u>(a-bc)/2  at  t=t^ 

*(t‘)  =' 

I 

La(ti)-2u]/b 

if  u<(a-bc/2  and  a<bh  at  t=t^ 

(33) 

Mt)  =  x(fw+)+h(t-tfw) 

for  all 

(34) 

u<(a-bc)/2  and  a>bh 

for  all  te(t( J) ,tfw) 

(35) 
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-W</J  [u-(a-bp) ]dt<0  for  all  te(tfw,t.) 

tfw 


(36) 


-W =  /*  [u-(a-bp) ]dt 
tfw 


(37) 


If  all  above  constraints  have  one  or  more  solutions,  then  we  save 
them  as  successor  vertices  of  v.  and  go  to  Case  5.2.  Otherwise, 

sJ 

terminate  this  branch  and  go  to  Case  5.2. 


Case  5.2  Let  ti  be  the  entry  time  to  I(t)=W.  This  case  is 
similar  to  Case  5.1,  except  that  (33)  and  (37)  are  replaced  by 
(38)  and  (39),  respectively. 

x.(t“)  =  [a(t.  )-2u/b  if  u<(a-bc)/2  and  a>bh  at  t=t  (38) 
i  i  -  l 


0=7.*  [u-(a-bp)dt 
^fw 


(39) 


If  there  exists  some  solutions  to  (t^,^)  in  Case  5.1  cr  Case 
5.2,  then  we  keep  them  as  successor  vertices  of  Vj  and  go  to  Step 
2.  Otherwise,  v.  is  fathomed  and  to  to  Step  4. 

v 

Step  6  Assume  that  t is  the  exit  time  from  l(t)  =  0,  and  ^  is 
the  next  entry  time  to  1  =  0  cr  I=W,  with  tj.Q<t^. 


Case  6 . 1  Suppose  that  ti  is  the  entry  time  to  1=0.  Solving 


the  following  constraints,  we  may  get  the  values  of  tfg  and  tj. 


,c 


X(tfo  )  s 


if  u>(a-bc)/2  at  t=t 


fO 


(40) 


'•[a(tf>0)-2u]/b  if  u<(a-bc)/2  and  a<bh  at  t=tfQ 


i(t.“)  =■( 


if  u>(a-bc)/2  at  t=ti 


[a(t^)-2u]/b  if  u(a-bc)/2  and  a<bh  at  t=t^ 


(41) 


i(t)=X(tf0+)+h(t-tf0+)  for  all  te(tf0,ti) 


(42) 


0</fcl  [u-( a-bp) ]dt<W  fcr  all  te(t~n,t.) 
t  fO  i u  i 


(43) 


Or/.1  [u-(a-bp) ]dt  (44) 

*fO 

u>(a-bc)/2  cr  [u<(a-bc)/2  and  a<bh]  fcr  all  te ( t^ ^ ,t^ q)  (45) 


We  keep  all  solutions  as  successors  of  v  ,  if  any,  and  go  to  Case 

6.2. 


Case  6.2  Suppose  that  is  the  entry  time  to  I=W.  Again, 
this  case  is  similar  to  Case  6.1,  except  that  we  replace  (41)  and 
(44)  by  (46)  and  (47),  respectively. 

x(ti')s[a(ti)-2u]/b  if  u< ( a-bc)  /2  and  a>bh  at  ts^  (46) 

W =/  [u-(a-bp)]dt 

fO 


(47) 
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If  there  exist  any  solutions  to  ^ase  ^.1  or  Case  6.2 

then  save  them  and  go  to  Step  3*  Otherwise,  check  whether  (45)  is 
satisfied  by  setting  tf0=T.  If  yes,  let  l(t)=0  on  [t(j),T]  be  a 
feasible  solution  and  go  tc  Step  2.  If  not,  Vj  is  fathomed  and  go 
to  Step  4. 

Step  7  (Termination)  J  is  optimal. 

To  illustrate  this  algorithm,  we  will  solve  a  simple 
numerical  example. 

An  Example 

Suppose  that  T=10,  a(  t) =30+1 Ot-t^ ,  b=1,  h= 1 ,  c=10,  u=10, 

W=125/12  and  IQs102/12.  We  then  have  that  u<(a-be)/2  for 
t e ( 0 , 1 0 ,  a>bh  for  t<4.5  and  a<bh  for  t>4.5. 


Step  0  J=0 

Step  1  Solving  Case  1.1,  we  have 
x(t1")=10+10ti-t^,  tt<4 

x(t)=10+9ti+t-t^,  for  all  t<ti 

t, 

W-I0=4=/01[u-(a-bp) 3dt 
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which  gives  a  solution  t^=1,  x(t)=l8+t  for  t<1.  Let  the  vertex 

corresponding  this  case  is  v^  Next,  we  find  that  there  are  no 
feasible  solutions  to  Case  1.2,  and  go  to  Step  2. 

Step  2  t^sIjMO  and  go  to  Step  3. 

Step  3  J[0,1]+J  >0  so  that  v1  is  alive.  Go  to  Step  4. 

Step  4  There  is  a  unique  live  vertex  v^  with  1(1 )=W.  Go  to  Step 
5. 

Step  5  Solving  (32)— (37)  simultaneously,  we  obtain  tfw=2  and 
^=7.  Let  v2  be  the  vertex  corresponding  to  this  case.  Again,  we 
find  that  there  are  no  feasible  solutions  to  Case  5.2.  Go  to  Step 
2. 


Step  2  t^=7/10  and  go  to  Step  3. 

Step  3  J[0,7]+J7>0  so  that  v2  is  still  alive.  Go  to  Step  4. 

2 

Step  4  There  exists  only  one  live  vertex  v2,  and  l(t^)=0.  Go  to 
Step  6. 

Step  6  There  are  no  feasible  solutions  to  both  Case  6.1  and  6.2 
since  u<(a-bc)/2  and  a<bh  for  all  t>7.  We  know  that  I(t)  =  0  on 
[7,10]  is  a  feasible  solution,  let  v^  be  the  corresponding  vertex 
to  this,  and  go  to  Step  2. 
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Step  2  t[3)=10. 

J[0, 10]=/g°[p(a-bp)-hl-cu]dt, 

where 


1 

.18+t 

if 

0<t<  1 

10+10t-t2 

if 

1  <t<2 

♦<t)  = 

24+t 

if 

2<t<7 

M0+10t-t2 

if 

7< t< 1 0 

u(t)=10  for  all  0<t<10 

■24+5.5t-0.5t2 

if 

0<t<  "I 

20+ 1 Ot-t2 

if 

1  <t<2 

p(t)  = 

27+5.5t-0.5t2 

if 

2<t<7 

■  20+1 Ot-t2 

if 

7<t<10 

,102/12+4t-9t2/4+t3/6 

if 

o^to 

125/12 

if 

1<t  <2 

I  ( t)  = 

49/12+7t-9t2/4+t3/6 

if 

2<t<7 

1 

>0 

if 

7<t<10 

Go  to  Step  3. 


Step  3  J= J [0 , 1 0]  so  that  is  fathomed.  Go  to  Step  4. 

Step  4  There  are  no  live  vertices.  Go  to  Step  7. 


-30- 


Step  7  J=J[ 0,10]  is  the  optimal  solution.  This  shows  that  t =2  is 
a  maximal  strong  planning  horizon  and  t=7  is  a  minimal  strong 
forecast  horizon,  see  Figure  5. 


3 •  Strictly  Convex  Non-negative  Production  Cost  Case 

In  this  section,  we  extend  Pekelman's  model  [7]  by 
considering  the  production  and  warehouse  capacities.  We  will  see 
that  the  optimal  inventory  trajectories  and  the  strong  planning 
and  strong  forecast  horizons  for  the  linear  production  cost  and 
strictly  convex  production  cost  cases  are  quite  similar. 

Suppose  that  monopolist  faces  a  strictly  convex  non-negative 
production  cost  f(u)  instead  of  a  linear  production  cost  cu.  Then 
the  monopolist  will  attempt  to  maximize  the  mathematical 
expression 


^2=^0 a"kp) -hl-f( u)  > d t 


(48) 
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1 


subject  to  (2)  and  (3).  Here,  we  assume  that  f(u)  is  twice 
differentiable,  and  f'(0)<a/b  for  all  t. 

The  necessary  conditions  for  an  optimal  solution  to  (48)  are 
similar  to  those  of  (1)  except  that  Equation  (6)  is  replaced  by 
the  following  equation 


if  x+o1-p2<f' (0) 
if  f,(0)<x+p1-p2<f,(u) 
if  f* (u)<x+Pl-p2 

where 

g=(f,)‘1. 


* 

u  = 


,0 

g(  x  +  p- 
u 


(49) 


(50) 


3 • 1  Optimal  Policies 

Again,  there  are  only  three  different  possible  cases  for  the 
values  of  I(t)  in  this  problem  but  the  situation  of  optimal 
policies  is  slightly  different  to  the  linear  production  case. 

Case  1  0<I(t)<W.  This  implies  Pi=p2=0.  We  have  three 

regimes  la,  1b,  and  1c  in  this  case. 

Regime  la.  x(t)<f'(0).  We  yield 

u*=0,  p*=max{0;(a/b+x)/2} ,  i<0  (51) 

Regime  1b.  f '  (OX x(t)<g*  (u) .  We  obtain 
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u»=g(x)<0,  p^sminta/b;  (a/b+x)/2}  (52) 

Regime  1c.  f'(u)<x(t).  We  get 

u*=u,  p*smin{a/b;  (a/b+x)/2>  (53) 

Case  2  l(t)=0  and  l(t)=0.  This  implies  p2=0,  u  =a-bp*,  and 

g(x  +  p1)-Ca-b(x  +  p1)  ]/2  if  x+p^mintf '  (u) ;  a/b} 

t  ■  '  u-[a-b(x+p1)]/2  if  fXuXx+pXa/b  (54) 

>g(x+p.|)>0  if  a/b< x  +  p i < f '  (u )  . 

Thus,  we  have  two  regimes  2a  and  2b  in  this  case. 

Regime  2a.  i|»g  =  x+p^-p2<niin{f  ’  (u) ;  a/b}.  We  have 

u*  =  gUg)  and  p*=(a/b+4»Q)/2  (55) 

and  satisfies  the  following  equation  (56)  and  constraint  (57): 

g(*0)-a/2+b*0/2=0  (56) 

and  ^gsx+p^h,  or 

*0sa/[b+2g’(*0)]<h  (57) 

Regime  2b.  i»ga+p1-P2s(a-2u)/b>f '  (u)  and  a<bh. 
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»  .  * 

We  get  u  =u,  p  =  (a-u)/b,  u ^  =  n 1  =  n2=  0 »  and 
U£= ( a-2u-bc) /b>0 

Case  3  I(t)sW  and  i(t)  =  0.  This  implies  p^O,  u*  =  a-bp* 
also  have  two  regimes  3a  and  3b  in  tihs  case. 

Regime  3a.  it/w=\  +  p1-p2<niin  {f '  (u) ;  a/b>  and  a>bh . 

We  obtain  uSgt^),  p*  =  (a/b+>|i  )/2, 

and  i|»  satisfies  the  following  constraints: 
w 

g(4>w)-a/2+b*w/2=0 

and 

l>w=a/[b+2g' (*w) ]>h 

Regime  3b.  h»w=  x-«- p  1  - p2=  (  a-2u)/b>f '  (u)  and  a>bh . 

We  yield  u*=u ,  p*= (a-u)/b ,  u1=n1=n2=0,  and 
u2= ( a-2u-bc)/b>0 

Note  that  Case  3  can  occur  only  if  a>bh.  That  is, 


(58) 

.  We 

(59) 

(60) 

(61) 

(62) 

full 


warehouse  can  happen  only  if  the  rate  of  price  increase  is  not 
less  than  the  inventory  cost. 


3.2  Theoretical  Results  and  Planning  Horizon  Theorem 

Using  argument  similar  to  that  in  Section  2  it  is  clear  that 
Theorem  1  is  still  the  in  the  strictly  convex  production  cost  case 
but  Theorem  4  is  not  correct  in  that  case.  We  also  need  a  slight 
modification  of  the  arguments  to  get  results  similar  to  those  in 
Theorem2  2  and  3. 

Theorem  8.  If  Ig>T0  or  a/b<h,  then  i<0  for  all  t. 

Proof .  If  Iq>Iq  ^en  t*le  Pro°f  i-s  trivial.  Suppose  that 

—  • 

Iq<I0  and  there  exists  an  open  interval  (t,  t+e)  such  that  l(t)>0 

on  (t,  t+e).  Then  this  can  occur  only  in  Regimes  '1b  and  1c. 

Therefore,  we  have 

if  f ' (0)<x<min{a/b;  f'(u)} 
if  a/b<  x<f ' (u_ 
if  f ’ (u)<\<a/b 
if  x>max{a/b;  f'(u)} 

if  f ’ (0)<x<min{a/b;  f'(u)} 

if  f ' (u) <x<a/b 

Since  f(u)  is  strictly  convex  and  non-negative,  we  know  f'  and 
g-(f')"’1  are  increasing.  Hence,  gr  must  be  positive.  So, 
I( t»e)>0  which  also  implies  I(T)>0.  This  contradicts  that  fact 
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Theorem  9«  If  a/[b+2g'3>h  for  all  t,  then  l(t)>0  almost 
everywhere  on  [0,T]. 

Proof.  If  not,  l(t)=0  on  some  open  interval  (t,  x+e).  Since 
g'  is  positive,  we  have  a/b>a/[b+2g' ]>h  and  Regimes  2a  and  2b 
therefore  cannot  happen.  This  leads  to  a  contradiction. 

By  using  a  proof  essentially  identical  to  that  of  Theorem  7 , 
we  can  show  that  the  same  Planning  Horizon  Theorem  which  holds  for 
the  linear  cost  case  (Theorem  7)  also  holds  for  the  convex  cost 
case.  We  restate  the  theorem  for  completeness. 

Theorem  10.  (Strong  Planning  Horizon  Theorem) 

f  ** 

Let  t  be  the  exit  time  from  I(t)=W  and  t  be  the  next  entry 

time  to  l(t)=0.  Then  tf  is  a  maximal  strong  planning  horizon  and 

s» 

t  is  a  minimal  strong  forecast  horizon.  The  theorem  also  holds 
when  0  and  W  are  interchanged. 

4 .  Conclusion 

In  this  paper  we  have  studied  a  general 
price-production-inventory  model  with  linear  and  nonlinear 
production  costs.  We  characterized  the  optimal  trajectories  and 
showed  that  there  could  be  strong  planning  and  forecast  horizons. 
We  presented  a  forward  branch  and  bound  algorithm  which  identifies 
strong  planning  and  forecast  horizons,  and  uses  them  to  decompose 
the  problem  into  a  set  of  smaller  problems.  The  algorithm  is 
illustrated  by  means  of  a  simple  example. 
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It  would  be  possible  to  extend  the  results  of  this  paper  in 
several  different  ways.  For  instance,  more  general  demand 
functions  could  be  considered.  In  a  subsequent  paper  we  intend  to 
study  the  model  when  backlogging  is  permitted.  Also,  if  we  change 
the  linear  inventory  cost  into  a  strictly  increasing  non-negative 
inventory  holding  cost,  then  our  strong  planning  horizon  theorem 
is  still  true.  The  other  theorems  in  this  paper,  with  suitable 
modifications,  would  also  be  true. 
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separately,  and  their  solutions  put  together  to  form  a  complete  solution  of  the 
problem.  Finally,  we  derive  a  forward  branch  and  bound  algorithm  to  solve  the 
problem.  The  algorithm  Is  Illustrated  with  a  simple  example. 


